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We analyze electroweak baryogenesis during a two-step electroweak symmetry breaking transition, 
wherein the baryon asymmetry is generated during the first step and preserved during the second. 
Focusing on the dynamics of CP-violation required for asymmetry generation, we discuss general 
considerations for successful two-step baryogenesis. Using a concrete model realization, we illustrate 
in detail the viability of this scenario and the implications for present and future electric dipole 
moment (EDM) searches. We find that CP-violation associated with a partially excluded sector 
may yield the observed baryon asymmetry while evading present and future EDM constraints. 


I. INTRODUCTION 

The origin of the cosmic matter-antimatter asymme¬ 
try remains one of the outstanding mysteries at the in¬ 
terface of particle and nuclear physics with cosmology. 
The asymmetry is typically characterized by the baryon- 
to-photon ratio 

Yb = Pb/s = (8.59 ±0.11) X 10““ (1) 

where pB and s are the baryon number and entropy 
densities, respectively, and where the value has been 
obtained using data from Planck [T], WMAP [2] and 
large scale structure measurements. Assuming a matter- 
antimatter symmetric Universe at the end of the infla¬ 
tionary epoch, a non-vanishing Yb can be generated if 
the microphysics of the early Universe satisfies the well- 
known three “Sakharov conditions” [3]: (1) baryon num¬ 
ber violation; (2) C and CP symmetry violation; (3) de¬ 
parture from thermal equilibrium (or violation of CPT 
invariance). 

A variety of scenarios that satisfy these criteria have 
been proposed, each corresponding to a different epoch 
in cosmic history. One of the most widely consid¬ 
ered is electroweak baryogenesis (EWBG), wherein the 
baryon asymmetry. Ye, was generated during the era of 
electroweak symmetry-breaking (EWSB) that occurred 
roughly ten picoseconds after the Big Bang. Success¬ 
ful EWBG requires a strong first order electroweak 
phase transition (EWPT) that proceeds via bubble nucle- 
ation. CP-violating asymmetries generated at the bubble 
walls diffuse ahead of the expanding bubbles, catalyz¬ 
ing baryon number creation through electroweak (EW) 
sphalerons. The expanding bubbles capture the non¬ 
vanishing baryon number, which is preserved in the bub¬ 
ble interiors if the transition is sufficiently strong so as to 
quench the electroweak sphalerons. For a recent review 
of this scenario, see Ref. [4] and the references therein. 

In principle, the Standard Model (SM) contains all the 
elements needed to satisfy the Sakharov criteria in the 
context of EWBG. In practice, the value of the Higgs 
boson mass is too large to accommodate a first order 
EWPT; in a SM Universe, the EWSB transition is of a 


cross-over type[5]-[8]. Even if the Higgs boson had been 
sufficiently light, the GP-violating asymmetries associ¬ 
ated with GP-violating phase of the Cabibbo-Kobayashi- 
Maskawa matrix are too small to lead to the observed 
value of Yr jOUTT]. Thus, successful EWBG requires both 
new degrees of freedom to induce a first order EWPT and 
new GP-violating interactions to generate sufficiently siz¬ 
able asymmetries during the transition. 

In this study, we consider EWBG in light of the pos¬ 
sibility that EWSB proceeded in multiple steps, rather 
than in a single transition from an electroweak-symmetric 
to an EWSB vacuum as has been conventionally as¬ 
sumed. We further consider the possibility that Yb was 
generated during a transition to a EWSB-vacuum that is 
not the present day vacuum (“Higgs phase”), but rather 
one that was the lowest energy state for a period prior to 
the final transition to the Higgs phase of the SM. For sim¬ 
plicity, we concentrate on a two-step scenario, though the 
general features could generalize to patterns of EWSB 
that entail additional intermediate phases. 

In earlier work, we demonstrated the viability of the 
two-step EWSB scenario induced by the presence of an 
electroweak triplet Y ^ (1,3,0) [12], in which the tran¬ 
sition to the first EWSB-vacuum involving only a non¬ 
vanishing neutral triplet vacuum expectation value (vev) 
is strongly first order. Some general conditions needed 
for a successful two-step EWPT in extensions of the SM 
Higgs sector were subsequently studied in Ref. [13]. Here, 
we focus on the generation of GP-violating asymmetries 
during the transition to the penultimate EWSB-vacuum. 
Starting with Y-extended SM, we show that generation 
of these asymmetries requires additional field content. 
After discussing general considerations, we focus on a 
concrete example that provides a proof-in-principle of vi¬ 
ability of the general paradigm. 

It is interesting to ask about the experimental signa¬ 
tures of the multi-step scenario. Requiring that the final 
transition to the SM Higgs phase occurs at sufficiently 
low temperatures as to avoid baryon number erasure 
through re-excited EW sphalerons implies that at least a 
subset of the mass parameters in the Lagrangian are not 
too different from the EW scale. Initial studies of the con- 
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sequences for collider phenomenology in the S-extended 
SM are discussed in Refs. [T^llTi]. In general, the intro¬ 
duction of new CP-violating interactions must contend 
with severe constraints from searches for the permanent 
electric dipole moments (EDMs) of atoms, molecules, and 
the neutron (see, e.g. Refs. [TSHIZ]) as well as possi¬ 
bly probes of CP-violation (CPV) in the heavy flavor 
sector for some scenarios p!8l - [2Q) . It is reasonable to ex¬ 
pect that any new CPV as needed for successful EWBG 
will be testable with the next generation EDM searches 
or heavy flavor studies. In what follows, we show that 
this expectation may not be borne out for the multi-step 
scenario, as the CP-violating interactions may involve 
a partially secluded sector whose impacts on low-energy 
CP-violating observables are highly suppressed. In short, 
multi-step electroweak baryogenesis may open a new win¬ 
dow for generation of Yb at the weak scale, one that is 
relatively immune to experimental constraints on CPV 
in the near term. 


Before proceeding, we note that others have consid¬ 
ered baryogenesis scenarios going beyond the conven¬ 
tional paradigm of a one-step EWSB. Ref. [21] consid¬ 
ered a two-step phase transition (2SPT) scenario using 
a Two Higgs Doublet Model (2HDM), wherein the first 
step is a second order (or cross over) PT, while the second 
step is strongly first order. However, because electroweak 
symmetry is broken during the first step, the H + L vi¬ 
olating processes are suppressed, and the CPV asymme¬ 
tries generated during the second step cannot be effi¬ 
ciently transferred into the baryon asymmetry [22]. The 
authors of Ref. [23] considered an extension of the SM 
with a complex scalar singlet [24], in which the universe 
first undergoes a transition to an EW-symmetric phase 
with a non-vanishing singlet vev, followed by a transi¬ 
tion to the Higgs phase. CPV-asymmetries are induced 
during the second step by a non-renormalizable Higgs- 
singlet-top quark interaction. Our scenario differs quali¬ 
tatively from these earlier studies, since Yb is generated 
during the first EWSB transition and is preserved dur¬ 
ing the subsequent transition to the Higgs phase. In this 
initial study, we also focus solely on CPV in the scalar 
sector involving only renormalizable operators. 


Our discussion of two-step electroweak baryogenesis is 
organized as follows. In section |n| we outline general 
considerations for baryogenesis in this scenario. In sec¬ 
tion HI we define the details of our model with the par¬ 
ticle content, interactions, relevant Eeynman rules and 
the conditions for EWSB. Section [IV| gives the frame¬ 
work for implementing constraints from two relevant ob¬ 
servables: EDMs and Higgs to diphoton decay rate. In 
section |V| we present the set of transport equations that 
describe the dynamics of particle-antiparticle asymmetry 
generation during the first step of the two-step transi¬ 
tion. Einally we present our numerical results in section 
EH and conclude in section Em Technical details associ¬ 
ated with solving the transport equations appear in Ap¬ 
pendices. A reader interested primarily in the general 
framework, specific model realization, and primary re¬ 


sults may wish to concentrate on section [Til th e first part 
of section m and Eigs. mu of section EH that show 
the sensitivity of present and future EDM searches to re¬ 
gions of the model parameter space consistent with the 
observed baryon asymmetry. 


II. GENERAL CONSIDERATIONS 

In what follows, we will adopt a specific model real¬ 
ization to illustrate the viability of two-step electroweak 
baryogenesis. Before doing so, we provide some general 
considerations that should guide the choice of a model, 
concentrating here on the ingredients needed to gener¬ 
ate a baryon asymmetry during the first step. We note 
that generation and preservation of the baryon asymme¬ 
try during the first step requires that it be a strong first 
order EWSB transition, while the preservation of this 
asymmetry during the second step to the SM Higgs phase 
requires that (a) the temperature of the latter transition 
is sufficiently low as to avoid re-exciting the electroweak 
sphalerons and (b) the entropy released during the second 
transition be sufficiently small so as to avoid over dilution 
of the asymmetry generated during the first step. The 
possibility of satisfying these requirements was demon¬ 
strated in our previous work on the E-extended SM [12] . 

In order to produce a non-vanishing Yb during the 
first transition, one requires generation of CP-violating 
asymmetries that ultimately yield a non-vanishing num¬ 
ber density of left-handed (LH) fermions, ul- The latter 
biases the electroweak sphalerons in the unbroken phase 
ahead of the advancing bubble walls whose interiors con¬ 
tain the first broken phase. We consider two sectors for 
this purpose: (a) the SM and (b) a new sector that con¬ 
tains the fields responsible for the first EWSB transition 
- generically denoted (j)j - plus additional fields that in¬ 
teract with these fields and that may be partially or com¬ 
pletely secluded from the SM. The following possibilities 
then emerge: 

• The new sector contains additional LH fermions 
that contribute to the B+L anomaly. CP-violating 
interactions of these fermions with the <pj lead to a 
non-vanishing 

• A CP-violating asymmetry is generated for one or 
more of the new sector scalar fields (j)j and is sub¬ 
sequently transferred to the LH fermions of the SM 
through interactions between the two sectors. 

• A CP-violating asymmetry involving SM fields is 
generated through their interactions with the (j)j 
yielding a non-vanishing ul either directly or indi¬ 
rectly via SM Yukawa and gauge interactions. 

In the remainder of this initial study, we concentrate 
on the third possibility, as it presents the greatest poten¬ 
tial for experimental accessibility; we defer a considera¬ 
tion of the other possibilities for future work. We also 
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focus on the case where the initial CP-violating asymme¬ 
try is generated entirely in the scalar sector and trans¬ 
ferred to the SM LH fermions via Yukawa interactions. 
In this instance, one requires at least two distinct scalar 
fields that mix through their interactions with the bubble 
walls of the first EWSB transition. For concreteness, we 
will utilize a 2HDM for the Higgs sector, where the dou¬ 
blet fields mix during the first EWSB transition. During 
this first step, however, neither of the doublets obtain a 
vacuum expectation value (vev). Consequently, we may 
treat both the neutral and charged components of the 
doublets as complex scalars with masses determined by 
the finite-temperature potential and the space-time vary¬ 
ing vevs of the fields driving the first EWSB transition. 
This treatment differs from what is appropriate when the 
neutral components of the doublets obtain vevs, leading 
to one combination of the CP-odd neutral scalar that is 
eaten by the Z-boson to become its longitudinal compo¬ 
nent and the other combination that is a physical CP-odd 
scalar. The latter framework, wherein the neutral CP- 
even and CP-odd scalars are treated as distinct degrees of 
freedom, applies to the second EWSB vacuum, or Higgs 
phase. 

For the first EWSB step, then, we consider two com¬ 
plex scalars hi and ^2 whose mass-squared matrix has 
the form 




m?i mfs A 
^12 ^22 / 


( 2 ) 


S,S' 



FIG. 1: Illustration of the two-step PT that we study in this 
paper. We focus on the BAU generated during the first step. 


doublet vevs remain zero during the first step, possible 
forms of the CPV interactions that satisfy these consid¬ 
erations include 

(a) CPV-asymmetries generated entirely in the new 
scalar sector: 

(j)\(t) 2 (l>j<l>j+h.c. , j = l,2 (3) 

(b) CPV-asymmetries generated in the SM (doublet) 
scalar sector: 

+ ,etc. (4) 


As discussed in Refs. [25l [26], generation of non¬ 
vanishing hi and ^2 asymmetries requires that 771^2 
contain a spacetime-dependent complex phase, 0{x) = 
Arg ( 77712 )- The latter arises from the interactions of 
the hi ^2 with the expanding bubble wall that effectively 
provides a spacetime-dependent background field via the 
space-time variation of the vevs of the fields (j)j involved 
in the first EWSB transition. 

Additional comments regarding the origins of 0{x) are 
in order. In principle, this phase could arise entirely 
spontaneously due to the vevs of the <pj without any de¬ 
pendence on an explicit CP-violating phase in the La- 
grangian of the theory. In practice, the net asymme¬ 
try produced by a coalescing ensemble of bubbles admit¬ 
ting only spontaneously-generated CP-violating phases 
will be zero [27]. For every bubble having a spontaneous 
phase of a given magnitude and sign, there will always 
be a partner bubble somewhere in the ensemble with a 
phase of the same magnitude but opposite sign. The con¬ 
tributions of the two bubbles will then cancel after the 
bubbles coalesce and the transition completes. Breaking 
the energy degeneracy between would-be partner bubbles 
requires the presence of an explicit CP-violating phase in 
the theory [2^. 

Under these conditions, 0{x) will be space-time depen¬ 
dent only when two distinct fields with differing non¬ 
constant space-time profiles contribute to 777^25 even in 
the presence of an explicit phase in the scalar potential. 
For the two-step transition of interest here, wherein the 


Note that the coefficients of these operators may 
be complex, yielding the requisite explicit CP-violating 
phase. Note also that only one of the (pj need carry 
SM electroweak charges, as required for quenching of the 
electroweak sphalerons during the first EWSB transition. 
The remaining new sector scalar fields may be pure gauge 
singlets or charged under other symmetries that do not 
contain the SM electroweak symmetries as a subgroup. 


III. THE MODEL 

With the foregoing considerations in mind, we illus¬ 
trate the viability of two-step EWBG with a concrete 
model. As indicated above, for the SM sector we extend 
the theory to a 2HDM in order to allow for mixing be¬ 
tween two SM-sector states. For the new sector we choose 
(pi to be the real triplet E ^ (1, 3, 0) and take the second 
field p 2 to be a real singlet S. In principle, we could have 
chosen p 2 to be a non-singlet with respect to SM symme¬ 
tries. Our rationale for choosing a singlet is that (a) the 
dynamics of EWSB during the first step are relatively 
simple; (b) it is relatively straightforward to write down 
a potential for S wherein the singlet-vev is non-zero and 
varying with temperature during the first transition; (c) 
the presence of the singlet illustrates the phenomenolog¬ 
ical features associated with a sector that interacts with 
the SM sector only through Higgs portal interactions of 
the type in Eq. With these choices, we take the vevs 
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of the singlet and neutral component of the triplet to be, 
Vs and respectively. The corresponding potential is 


.2 / . .X 


I + ^45 ^4 


= -^(s.s) + ^(s.s) +^s^ + 

i a2s hIH2 (s-s) + ^a2sHlH2S^ + h.c. 


+ 


+aissS-S 5 + -a2ssS-S + V{Hi,H2). (5) 


Note that the quartic couplings a 2 E and a 2 s are in gen¬ 
eral complex. The real couplings and a 2 ES' have 

little bearing on the dynamics of CPV, whereas they will 
play a role in the EWPT. For purposes of keeping the 
number of parameters manageable for this initial study, 
we will set them to zero. 

For the 2HDM potential we consider the Z 2 -symmetric 
limit under which Hi and H 2 are oppositely charged, 
thereby alleviates the possibility of potentially dangerous 
tree-level flavor changing neutral currents (FCNCs): 


V{H^,H2) = ^[HlH^y - 

+A3 [HIH 2 ) + A4 



[HIH 2 ) [hIh,) 


(6) 


Note that the Z 2 -breaking operators • F and 

hIH 2 S^ generate divergent contributions to the operator 
H\H 2 at one-loop order, implying the need for a counter 
term 772 ^ 2 ^ 1^2 + h-c.. We retain freedom to choose the 
finite part of mi 2 such that its sum with the finite parts 
of the one-loop graphs is sufficiently small as to satisfy 
experimental FCNC bounds. As we discuss in Section 
[VH there will also be finite-temperature contributions to 
H\H 2 associated with i7][i72F-F and H\H 2 S^ that can¬ 
not in general be eliminated and that will lead to mixing 
between the two doublets in the early universe. 

The pattern of two-step EWSB arising from this po¬ 
tential is illustrated in Fig.[^ At high temperatures, the 
universe starts in the completely symmetric phase with 
all vevs set to zero [58]. At a critical temperature T^, 
the vacuum with non-vanishing 7 ;^ and vs becomes de¬ 
generate in energy with the symmetric phase; just below 
Ts, the first order transition to the “F phase” proceeds 
through bubble nucleation. This transition is denoted by 
“Step 1” in Fig. 12 At a lower temperature Tf/, the vac¬ 
uum with non-vanishing vi and V 2 becomes degenerate 
with the F phase, and as the universe cools further, a 
transition to the Higgs phase occurs (“Step 2”). In the 
F-extended SM, it was observed that the second transi¬ 
tion is also typically first order. However, EWBG is not 


viable in this step because the sphalerons that are active 
before the first transition have been quenched in the F 
phase [ 59 ]. 

We will then investigate the possibility of EWBG dur¬ 
ing the first step. As a preamble, we first analyze the 
structure of the potential in greater detail. After EWSB 
occurs in the second step, the neutral components of each 
of the Higgs fields acquire a vev, and the fluctuations 
around this value can be characterized by charged ), 
CP even and CP odd (A^) fields respectively: 

( Hi \ 

= I 1 , where 7 = 1,2. (7) 

The potential E(i7i, 7 ^ 2 , F, S) has three complex cou¬ 
plings: A 5 , a 2 E 5 An overall phase in these couplings 
is unphysical and can be rotated away via a rephasing 
transformation on the complex scalar fields 

H2=e^'^^H'2, S = E', S = S'. {%) 


A global phase O 2 — 0i can be absorbed into the following 
redefinition of the couplings and vevs: 




«2E = 


a2S , 


{viv^y = 


a2s — e 


'Ci2S , 


(9) 


Without loss of generality, we assume vi = and 
V 2 = 1 772 1 e*^. The second equation in the first line leads to 
^ ^ 0i — O 2 rephasing transformation on the sponta¬ 
neously generated phase. The transformation in Eq. ^ 
and Eq. leaves the Lagrangian unchanged and, thus, 
phases in the couplings that can be eliminated with such 
redefinition are unphysical. The model then contains the 
following three physical phases: 


(5s = arg [ttjs ^ 1 ^ 2 ] , 

6s = arg[a2si'ii'2] > 


( 5 a 5 = arg Ag {viv^Y 


( 10 ) 


Due to the p parameter constraints, at the zero- 
temperature triplet vev must be small. In what follows, 
we set it to zero: (F) = 0. For simplicity, we will also 
take the zero temperature vev of S to vanish. Given the 
symmetry of the potential under S ^ —S and F ^ — F, 
both S and the neutral triplet may contribute to the dark 
matter relic density in this case. 


A. Minimizing the potential 

Stability of the vacuum state, after the transition to 
the second EWSB-vacuum, requires that the Lagrangian 
couplings and the vevs of the Higgs bosons satisfy the 
minimization conditions on the potential and positivity 
of all the masses in the spectrum. For E(i7i, 7 ^ 2 , F, S') 
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given in Eq. <§ and Eq. ^ we obtain the following min¬ 
imization conditions: 


mil = [Aicos^/3 + (As + A 4 + IA 5 I cos(5A5)sin^;d] 

+ cosfe + \a 2 s\vl cos(5^) tan/3, 

77222 = [^2 sin^ -h (As + A 4 -h IA 5 I COS 6 x ^) cos^ f3] 

+ (|<^ 2 e|'^I: COS(5e + \ci 2 s\vl cos 6 s) cot/3, 

sin/3cos/3 / /q 2 ix 1 • r 

0=-^- \^smp cosp V jAsIsmOAg 

sinc^E + \(i2s\vs 

0 = 7;e (^ 4 '^e — A^I + |<^ 2 E 1'^^ cos Je sin (3 cos /3) , 

0 = {b 2 s + + \ci 2 s\v‘^ cosSs sincos, 

( 11 ) 

where v = \/\vi\‘^ -h |'^ 2 p 5 tan^d = | 7 ; 2 |/| 7 ;i|. In this pa¬ 
per we concentrate on the case of no spontaneous CP 
violation^ i.e. ^ = 0 for simplicity. In this case the three 
physical phases are related to the complex couplings in 
the Lagrangian in the following way 


(3e = -arga 2 E, 6 s = -arga25, = -argA 5 .( 12 ) 


Note, that the phases 6 are manifestly rephasing in¬ 
variant, while the arguments of the complex couplings 
A 5 ,a 2 E 5 <^ 2 S' are not. The expressions above apply in the 
rephasing basis corresponding to the ^ = 0 choice that 
we have made. 

In the limit 7 ;e = = 0 we obtain minimization con¬ 

ditions that are identical to those of the 2HDM [29] 

772ii = [Aic^ + (A 3 + A 4 + Re A 5 ) , 

77222 = [^ 2^1 + (As + A 4 + Re As) c^] , 

0 = —Cj^Sj^ ImAs , (13) 


where = cos/3, s^ = sin/3 . The fourth and the 
fifth equations in Eq. (11), which correspond to equating 
to zero the partial derivative with respect to the triplet 
and singlet vevs correspondingly, become a trivial “zero 
equals to zero” equations. Note, however, that prior to 
the second step of the two-step transition, the latter two 
equations must be satisfied for non-vanishing and vs 
while 771 = 0 = 772 - 


B. Mass mixing 


do not mix and have a diagonal mass matrix 



ml+ 

0 


0 

to|+ 


(16) 


where the masses are related to the potential parameters 
via the following equations 


'^H+ = \ (“-^4 - Re A 5 ) v^, 

to|+ = -/x| + Re a 2 E c^sp . (17) 

The mass formula for the charged Higgs, ^ agrees with 
the corresponding Z 2 symmetric limit of Ref. m and is 
not modified by the presence of additional fields of the 
triplet and the singlet. 

The neutral scalar bosons mix and 

in the most general case (^ 7 A Q) there is one massless 
neutral Goldstone boson . In the case of our interest 
^ = 0 this state is the following combination of CP odd 
Higgs bosons A^^A^: 

= C/3H5 + 5/3H2 . ( 18 ) 


This equation is unchanged compared to the pure 2HDM 
Ref. |29| . Considering the mixing between the orthogonal 
state A^ = —spAiPcpA^ and three other neutral scalars: 
i ^2 5 ^^5 we obtain the following mixing matrix 

TAneutral 

/ Aic^ A345 0 ^ \ 

2 A345 SpCjB A2S^ 0 0 

^00 -ReAs 0 

\ 0 0 0 ^ / 

(19) 

where 772|o = ^ 1 + EO] 


A 345 — A 3 + A 4 -|- Re A 5 . 


( 20 ) 


The top left 3x3 block of the matrix is the 

same as in the scenario of pure 2HDM [29| for 772^2 = 0. 
Thus, there is no mixing between the triplet and the two 
Higgs doublets in our theory for both charged and neutral 
states. 

The singlet in our theory does not mix with any other 
neutral particles and its mass in terms of the parameters 
of the potential equals 


m\ = & 2 S + Re a 2 s C/ 3 S /3 . 


( 21 ) 


We start from the mixing among the charged particles 
i 7 ^,i 72 ^,E+. The massless Goldstone combination is 

G+ = C 0 H+ + S 0 H+ . (14) 

The remaining two orthogonal charged scalars 

4 ’l = -S/S-ffl" + C13H} , 

0+ = E+, (15) 


C. Relevant Feynman rules 

We concentrate on the type-H 2HDM, which is moti¬ 
vated by the minimal supersymmetric extensions of the 
SM. It has the following interaction Lagrangian between 
the Higgs bosons and the fermions 

C\i = —Yu Ql '^ 0-2 H 2 Ur — YdQl Hidn - 1 - h.c. ( 22 ) 









6 


TABLE I: Table of the parameters in the potential versus 
the phenomenological parameters. 


Parameters in the potential 

Phenomenological parameters 

Ai , A 2 , A 3 , A 4 , ReAs 

mh , 171^2 , Rea 2 s , Ima 2 s , Ms 

64 s , Rea 2 s , Ima 2 S' , b2s , b^s 

V, tan /3, a, Re a 2 S, Re a 2 s 

6^, 6s, ^4s , b^s , rri'E , ms 

5 ? nr/j,2 5 T7ih^ 


For the EDM constraint we will require the Feyn¬ 
man rules for the Yukawa interactions of the Higgs 
bosons with the fermions hiff, the tri-scalar interactions 
and the couplings with neutral gauge bosons 
yY+Y”. The relevant interaction Lagrangian 

reads 


r — J 


m/ 

V 


hi {cf,iff + Cf^i /ns/) + df /ns/ 


+ [S"*" {idfiTi ) — (i5^E~’’) S ] (e + <72 cw Z/j,) , 

(23) 

where the couplings c, c, A equal to 


^t,i — -^ 22 /<^/3: 

dt = - 1 , 


Cb,i = Ril/CjS, 

dd ~ I 5 


(24) 


Xi = {Rii S/3 + Ri2C^) Rea2E — Ri3^^CL2j: • 


In equations above the matrix R is defined 
to diagonalize the neutral bosons mass matrix 
-R^neutrai^^ = diag(m^^, J. In terms 

of matrix R the weak eigenstates are related to mass 
eigenstates via = {hi, h 2 , hs, h^) • R. 

Note that the fermions directly couple to the neu¬ 
tral Goldstone boson G^, while the scalar interaction 
Y+Y^G^ is absent at tree level. 


D. Phenomenological parameters 


phenomenological parameters via m 

Ai = 

A 2 = 

A3 = 

A 4 = 


rnl. 


ml^ 

tan^ [3 



y2 

1 

ml. 


ml. 

cot^ P 



v2 


ml. 

- 

ml. 

+ 2771^ 



v2 



- 

2 m^+ 


Re A 5 — — 


m 


hs 


2 2 2 
^11 = ^22 = 


(25) 


The matrix Rij that enters the Feynman rules in 


Eq. (24) in general is a function of three angles a,ab,o:c 

)■ 


R { ac , ab , a ) = R23(<^c)Ri3(<^5)^12 
In the case of interest for us we have 

/ C/3 S/3 0 0 


7r> 

2. 


R(0,0,^-^/2) = 


-sp C/3 0 0 
0 0 10 
Vo 0 0 1 


(26) 


Finally /i| is found from Eq. 0 , b 2 s is found from 
Eq. (21), and the imaginary parts of a 2 ^,a 2 s are found 


from their real parts and the angles S^,5s • 

The formulae described in this subsection allow us to 
recover all the parameters in the potential starting from 
phenomenological parameters in Table [T] . 


IV. OBSERVABLES 

In this section we review two observables that place 
constraints on our illustrative scenario: EDMs and the 
Higgs diphoton decay rate. The former constrains the 
interplay of masses and the CP-violating phase (5^. The 
latter is sensitive to the magnitude of Higgs portal cou¬ 
pling a 2 E and the triplet masses. 


We are interested in the zero-temperature relations 
among the physical parameters of the theory such as 
masses of neutral and charged scalars, and the param¬ 
eters in the potential. In Table |T] we list the set of pa¬ 
rameters of our potential and the phenomenological pa¬ 
rameters. 


In the 2HDM sector the CPV is absent because we take 
m ^2 = O 5 therefore ImAs = 0 from the minimization 
conditions in Eq. (13). The only two CPV phases in our 
theory are due to the triplet and singlet and are repre¬ 
sented by Sji, 6s- 

Motivated by present fits on to Higgs observables (see, 
e.g., [29] ), we assume the SM alignment limit a = 
P — IT 12, in which the couplings Ai,... A 5 are related to 


A. Electric Dipole Moments 

EDMs of non-degenerate systems are CP-violating ob¬ 
servables, and EDM searches provide constraints on BSM 
sources of CPV. Treating the SM as an effective field 
theory, new sources of CPV can be characterized by 
dimension -6 operators, which include elementary fermion 
EDMs, quark chromo-EDMs, the CPV 3-gluon operator, 
etc. (see review HE]). In general, a model-independent 
analysis of EDM limits would require us to consider a 
large set of these CPV operators. Here, we can focus 
on the electron and quark EDM operators, as these are 
the only experimentally relevant dimension -6 operators 
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where the dimensionless tensor charges have been 
evaluated in Ref. m and we use them below in section 

MM 

While the down EDM is proportional to tan p like the 
electron EDM, the up EDM behaves as cot (3 due to the 
different Yukawa structure. This means that the neutron 
and proton EDMs can potentially place a stronger bound 
on the low tan (3 region in comparison with the electron 
EDM bound. 


FIG. 2: The dominant EDM contribution comes from these 
Barr-Zee diagrams and their mirror graphs. 


that can be generated in our model at lowest non-trivial 
(2-loop) order [62]. 

We have assumed that our 2HDM sector is CP- 
conserving, and that the new sources of CPV involve ei¬ 
ther the triplet or the singlet. As a result, the leading 
contribution to SM fermion EDMs comes from graphs in¬ 
volving the charged scalar E+ shown in Eigurej^ These 
so-called Barr-Zee diagrams are familiar from EDM anal¬ 
yses in 2HDM, and the result for dimensionless EDM 
5f = —v^df/2mfe (as in [16]) from this diagram is anal¬ 
ogous to the loop result [30] : 


(i/jsv= 


2.,2 


Qfe^v 


2567r^m|+ 

,, _ 9zff 92 cwv^ 

^ ~ 2567r4m|+ 




(27) 




i=l 


where the vector coupling of the electron to the Z equals 
Qzee = 92 {sw “ l/4)/cvi/ and 

4=«^I+ML- (28) 


The loop functions are listed in Ap¬ 

pendix A. Analyzing expressions for in Eq. (24) in the 
limit ab = ac = m \2 = ImAs = 0, we see that only 
i = 3 contributes to the fermion EDM. Thus theoretical 
predictions for both electron and neutron EDMs depend 
only on the masses rrih ^, m^+ . The electron EDM which 
gives the strongest bound is proportional to the sin 
and tan p parameters: 5e ^ sinc^^ tan p. 

Up and down quark EDMs contribute to the neu¬ 
tron/proton EDM as 


dn/p — ^ Cn i-dn 




(29) 


q=u,d 


where are matrix elements of quark EDM operators 
between neutron/proton wavefunctions. They equal 


^n/p ..2 ’ 


(30) 


B. Higgs Diphoton Decay 


CMS [32] and ATLAS [33] measurements of the Higgs 
diphoton decay are consistent with SM and provide 
bounds a 2 E and ms- The observed signal strength for 
the Higgs diphoton decay channel, divided by the cor¬ 
responding SM value, is defined as g. Approximating 
the SM value of the diphoton decay signal strength with 
contributions from the top and VF—boson loops, using 
standard results from Ref. [34] we obtain the following 
theoretical prediction for the value of g 


theory _ 

/^2HDME ~ 


\^9t Q‘tPi/ 2 {Tt) + gw Fi{tw) + 5s -Fo('rs)| 

|3Q|Fi/2(Ti) + Fi(TH/)|' 


(31) 


where T( = 4m‘f/ml^,Tw = 
and [34] 


Fxl2ij) = -2t[ 1 + (1 - T)/(r)], 

Fi(t) = 2 + 3t[ 1 + (2 - r)/(r)], 

Fo{t) = t[1 - r/(r)]. (32) 


Eunction /(r) reads [34] 


f(.T) 


&rcs\v?{l/^/ t), if r > 1 

“I (injUr^ “^’^) 


(33) 


The coefficients gt^gw^dT. are 


gt = 1 , 

gw = 1, 

^ sin 2/3 Re a 2 E 



(34) 


The numerical constraints from the Higgs to diphoton 
data on the parameters tan/3,mE are presented in the 
section IVIBI below. 


V. BARYOGENESIS AND TWO-STEP PHASE 
TRANSITION 

We now analyze the dynamics of BAU generation dur¬ 
ing the first step of the 2SPT process. Eollowing early 














work, we first compute the net number density of LH 
fermions ul generated by CP-violating interactions at 
the bubble walls (the space-time varying scalar field 
vevs). As it diffuses ahead of the bubble wall, the non¬ 
vanishing riL then catalyzes creation of a non-vanishing 
baryon number density by the EW sphalerons. Since the 
EW sphaleron rate E^s is typically much slower than the 
rates for processes that govern ul generation, we treat 
the computations oi pB and ni, as separate steps. The 
latter entails deriving and solving a coupled set of trans¬ 
port equations of the form: 

dxPk = -Y.^A{l^k-l^i----) + sf , (35) 

A 

where and pk are the number density current and 
chemical potential for particle species “/c”, E^ are a set 
of particle number changing reaction rates that involve 
species k and other species relevant to the problem, and 

is a CP-violating source for species k. 

We derive these equations using the Schwinger- 
Keldysh closed time path (CTP) framework that is ap¬ 
propriate for out-of-equilibrium dynamics at finite tem¬ 
perature [35H4Q] . Eor a detailed review of the formal¬ 
ism in the context of the EWBG in SUSY see e.g. 
Refs, [m 02] . Eor brevity we provide the main results 
here without going into the details of the formalism. 
As this work provides an initial study of the viability 
of two-step EWBG, we employ several approximations 
and assumptions to make the computation reasonably 
tractable, deferring a more exhaustive treatment to fu¬ 
ture work: 

• We take the bubble walls to be planar, and treat 
all quantities only as functions of the co-moving co¬ 
ordinate z = X + where is the wall velocity 
and X is the position relative to the center of the 
wall. We illustrate the dependence of Yb on 
and give illustrative results for a value of that 
is within the range of values obtained from EWBG 
studies for other models. 

• We use bubble wall profiles, vj:{x) and vs{x), that 
have the typical form obtained in other studies of 
bubble walls, again showing the dependence of Yb 
on the parameters that characterize the profile and 
choosing typical values for purposes of illustration. 


The VIA assumes that the particle-antiparticle asym¬ 
metry generation is dominated by the region near the 
phase boundary, where the vevs are small compared to 
both T and the difference \m\i — 77122 While the ap¬ 
proximation is thought to provide a reasonable estimate 
of the magnitude of the GEV sources, it is associated 
with theoretical uncertainties. In particular, it neglects 
the impact of flavor oscillations, which become impor¬ 
tant in the region where the off-diagonal term in is 
comparable to 1777^;^ ~ ^ 22 1- present instance, the 

potential V{Hi,H 2 ) contains no tree-level contribution 
to 777^2 would give rise to flavor oscillations. Ela- 
vor oscillations will, nevertheless, be induced by a finite- 
temperature contribution as well as the non-vanishing 
vj: and vs- Thus, we will treat our results for parameter 
choices in this region with a healthy dose of salt. 

A more complete treatment that includes flavor os¬ 
cillations requires a resummation of the vevs to all or¬ 
ders, as well as a first principles treatment of the ther- 
malizing interactions between the scalar fields and the 
other particles in the finite-T plasma. Initial efforts to 
carry out the former have been performed in Refs. [13]- 
1^ . focusing on the Yb in the minimal supersymmetric 
Standard Model as generated by gaugino-Higgsino inter¬ 
actions with the bubble walls (for an extensive discus¬ 
sion and related references, see Ref. HD- The authors of 
Refs. 05] 06] found a significant reduction in the asym¬ 
metry compared to the result in the vev-insertion ap¬ 
proximation. However, as later pointed out in Ref. [26] . 
the computation of Refs. 05] 06] neglected the effects of 
diffusion ahead of the bubble wall and dropped the dom¬ 
inant GEV source for the flavor-diagonal particle number 
densities. Retaining the latter and including the effects of 
diffusion leads to a “resonant-enhancement” of Yb in the 
small \mii — 77722 1 regime that is consistent with what is 
observed in the VIA. Consequently, we conclude that for 
our present purpose of evaluating the viability of two-step 
EWBG the vev-insertion approximation provides a rea¬ 
sonable estimate of the magnitude of the baryon asym¬ 
metry one might anticipate in this scenario. We defer 
to future work an analysis using the more numerically 
intensive framework of Refs. [25l [26] that also requires 
modeling in detail the CE-conserving thermalizing inter¬ 
actions in the high-temperature plasma. 

A. Relevant interactions 


• We compute the GEV sources using the vev- 
insertion approximation (VIA), which amounts to 
expanding the mass-squared matrix to second or¬ 
der in the off-diagonal elements 777^2 (^)- The 
sources then depend on the interference of these 
elements at two different space-time points, viz, 
777^2 (^)^ 12 (^)*- will below, when the 

GEV phase in the 777^2 (^) vary with space-time, 
the interference 777^2 (^)^i 2 (^)* contains a non¬ 
vanishing GEV phase. 


With these comments in mind, we now provide the 
Lagrangian for interactions that are responsible for the 
BAU generation below. All scalar interactions are con¬ 
tained in the effective potential V{Hi, S), see 

equations Eq. ^ , Eq. ^ . The source terms are gener¬ 
ated by the interactions that have two Higgs fields Hi, H 2 
and two VEVs of the triplet and singlet fields v^,Vg. The 
corresponding Lagrangian is given by 


= 

^scalar 


(Q2S + Q2S Vs) 


HlH2 + h.c. 


(36) 


2 
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FIG. 3: VEV insertion Feynman graphs in the 2HDM with 
an addition of a triplet and a singlet. The VEV and the a 2 
coupling in the graphs correspond to the triplet and a singlet 
correspondingly. 


The CP-conserving interactions include both scatter¬ 
ing and particle number changing reactions. The former 
determine how effectively particle asymmetries generated 
at the bubble wall diffuse into the broken phase where 
the weak sphalerons are unsuppressed. We take these 
into account using the diffusion ansatz (see below). The 
leading particle number changing reactions are mediated 
by tri-scalar and Yukawa interactions. For the former the 
following interaction Lagrangian is needed 


jC 


Y 

scalar 


(a 2 s Vs + a 2 s vs S) 


HlH2 + h.c. 


(37) 


For the latter, the Yukawa Lagrangian in the type-II 
2HDM is given in Eq. (22). Keeping only terms pro¬ 
portional to the top Yukawa coupling we arrive at 

-^Lmion = yt {-^2 tRti + H 2 iRbs) , (38) 


where yt is related to the (3,3) component of the Yukawa 
matrix Yjj and the top quark mass in the following way 
yt = = V 2 mt/\ V 2 1 . The corresponding expres¬ 
sion for the bottom Yukawa coupling is ^5 = = 

\/2mi)/\vi\. Note that in the regime of moderate to 
large tan/3, wherein ^5 is enhanced over its SM value, 
explicit inclusion of bottom Yukawa interactions can be 
decisive [47]. 


B. The source term 

Under the VIA, the processes shown in Figure [^gen¬ 
erate the source term densities iLi,iL 2 - The complex 
couplings that enter the computation are a 2 E,<^ 2 S'- As 
a result the source term has two qualitatively different 
parts: the CP conserving and the CP violating. Apply¬ 
ing straightforwardly the CTP method [41], we obtain a 
CP conserving term 

{yH2 - y-Hi), 

3^2=-SW., (39) 

where the CP conserving relaxation rates T are 

Tjui = -^^^\a2EvU^) + a2svUx< ■ 

(40) 


xV{x) 


y{y) 


xv{x) 


y{y), 


\< H2 






Hi 


Hi 


\ 

Hi 


/ 

'H 2 


Ax) Hi v{y) 


\ 

H 2 


V H 2 


FIG. 4: First two graphs define the relaxation rate F/f. The 
last graph leads to a zero source term for the held. 


The role of the CP conserving term in the Boltzmann 
equations is that of the relaxation type. The CP violating 
source term equals 

Sf{x) = 

TT^ 

X [us(a;)uE(a:) - us(a;)uE(a;)] A. (41) 

The role of the CPV source term in the Boltzmann equa¬ 
tions is to generate the particle-antiparticle asymmetries. 
The thermal functions , A depend on the tempera¬ 
ture T and thermal masses of Higgs bosons Hi,H 2 and 
are given by 

_ f k‘^dk 1^^ / hB{e 2 ) T /iB(ef) _ hB{e 2 ) T hB{ei) \ 

J LO1LO2 2 \ €2 — Cl €2 + ei / 

_ r k^dk /ns(ei) — ns(e2) Tis(ei) + Tis(e2) + 1 \ 

J lj;^2 V (U - 62)2 + (e2 + ei)2 ; ’ 

(42) 

where — iF//., i = 1,2. The func¬ 

tions , A have a well known resonant enhancement 
when rriHi = hih 2 im, which we stress below in the 

Section EEl 


C. Particle Number Changing Rates 


The hrst two graphs shown in Fig. [^ evaluated us¬ 
ing the CTP formalism lead to terms in the trans¬ 
port equations for the currents of the densities Hi^H 2 '. 
3A Hi = Y Ahi - with i = l,2, where 


^ 2^2 


I tt2E I {vy (x) ; rriHi , 'rnH2 , ) 

1^25 {vs{x)]mHi,mH2,ms) .(43) 


The approximate sign indicates the assumption of slowly 
varying triplet and the singlet vevs in the first two graphs 
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section, we construct the following four densities: 

Q = nt^ 

h = n^+ + Ufjo + + Ufjo , (46) 

where T^Q^H^h correspond to the densities of right 
handed quarks, left handed quarks and two linear com¬ 
binations of Hi,H 2 densities, respectively. The resulting 
set of the transport equations are 


FIG. 5: Processes in the plasma induced by the Yukawa terms. 


of Figure]^ ~ • The function depends 

on the temperature and thermal masses of i^i, 7^2, S and 
S fields, it can be found in Ref. m and is not repeated 
here for brevity. The third graph in Figure equals to 
zero in the CTP formalism, which is manifestation of the 
fact that the field has a zero chemical potential. 

In Figure we show the additional interactions in the 
plasma that lead to Yukawa Ty and corresponding ad¬ 
ditional terms in the transport equations that arise from 
Yukawa interactions. These processes play important 
role as they transfer the particle-antiparticle asymmetry 
generated in the scalar sector to the left handed quark 
anti-quark asymmetry. The type II 2HDM has a Yukawa 
sector equivalent to MSSM, if one excludes the superpar¬ 
ticles from the latter Lagrangian. Thus the relaxation 
rate Ty for our scenario and the corresponding terms 
in the transport equations can be readily obtained from 
MSSM calculations [42]: 

12Arc ^ , n 1 on 

Ty = — — — If {mt^.mQ.rriH^) + 0-129— T , 

"" (44) 


= -Vy 


+r. 


= -Fy 


■ T 


f H hV 

ky 

kq 

\kH kh)_ 


^ ^ 9(Q + T) 

kq kr 


kg kr ' ' ’ 


-2r,, 


kq kr 


ks 

_h 

kn kh 
9(Q + T) 


ks 






-Fy 


'Q_ 

^ +1 

+-^y 

_kQ 


\kH kh) 


+ sf, 


= -ry 



1 -r ^ 

T ■ 

\kH kh j 


ky _ 


(47) 


The finite-temperature expressions for the coefficients ki 
can be found in Eq. (72) of Ref. [41] . For reference, we 
note that in the massless limit ky = kQ/2 = ks = 3 and 
kn = kh = 4:. Finite-temperature contributions to the 
thermal masses lead to modifications of these relation¬ 
ships. 


E. Approximate solution 


where the function Ip can be found in Ref. [42] and will 

2 

not be repeated here. The additional term 0.129|^T is 
included mi ss] as an estimate of the four-body con¬ 
tributions to the Fy, which plays role in some in the 
parameter space regions where, due to kinematic thresh¬ 
old effects, Xf vanishes. In particular, Xf vanishes for 
mH 2 <150 GeV, and for mH 2 = 200 GeV, the first term 
in the Eq. (44) ^ 3GeV. This to be compared to the 
value of the four-body term « 1.9 GeV. 

We use the standard form for the strong sphaleron 
rate 


Tss = l6Ka^ T, 


(45) 


with hi 1 and as is the strong coupling. 


D. Boltzmann equations 


Eor purposes of deriving intuition about the transport 
dynamics, it is helpful to proceed toward a solution to 
Eqs. (47) as far as possible analytically. To that end, 
it is useful to consider the limit Fy,Fss ^ oc, (see for 
example [41] and the references therein). In this limit 
the linear combination of densities that multiply Ty and 
Tss relax to zero very close to the bubble wall, yielding 
two conditions that allow us to eliminate two of the four 
densities. We note that one difference with [41] is that 
in our application we keep the h density in our coupled 
equations, thus when we use the approximation of large 
Fy, Tss we find 


^ _ H-^h 

kq kn 


0 , 


T 
ky 

2Q T ^ 9(g + T) _ ^ 

kq kr ks 


(48) 

(49) 


In order to analyze all the transport processes that 
arise from CTP processes considered in the above sub- 


Our four equations reduce to a set of two coupled Boltz¬ 
mann equations which we solve numerically. For that 
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reason we call this approximate solution as opposed to 
analytical solution which we would be able to achieve 
in the absence of the h density. Using the approximate 
formulas Eq. (48) and Eq. (49) we solve for T, Q. The 
solution reads 


T = CT{H^h), Q = CQ{H + h), 


where 


Ct = 


cq = 


kri^kB + 9^q) 
knikB + 9(^q + kr)] ’ 
kQ{kB — 9/ct) 
knlkB + 9 (/cq + kr)] 


(51) 


The remaining two densities i7, h satisfy a set of two cou¬ 
pled Boltzmann equations. In order to find these equa¬ 
tions we need to find linear combinations of equations in 
Eq. (47) that are free from either ry,rss . One possibil¬ 


ity is to choose 2 d^T^ + and — df^h^. 

The resulting equations should be solved with respect to 
H" ^ h" and as a result we obtain [63] 

Dh H" - {aiiH' + ai2h') - (fni/ + ri2h) + = 0, 

Dh h" - {a 2 iH' + a 22 h') - {T 21 H + fas/i) + ^2 = 0, 

(52) 

where we made the following definitions 


the vanishing of ri in this limit. A non-vanishing ub 
then arises from retaining the subleading terms in l/Tss, 
as contained in the r 2 -term in Eq. (84) of Ref. [41]. How¬ 
ever, keeping the finite temperature contributions to the 
quark and scalar field masses yields a non-vanishing re¬ 
sult at zeroth order in l/Tss- In Section |VI 


we compare 


(50) the solution of Eq. ( |5^ to the solution for the full set 


in Eq. (47) and we find for our benchmark scenarios that 
the approximate solution Yb is by ~ 10% larger than 
the result of the full solution. We studied the agreement 
between the approximate and full method when varying 
the parameters of the theory in the wide range. In some 
points in the parameter space the approximate Yb be¬ 
comes up to a factor of 2 larger than the full numerical 
solution, which cannot be explained by the subleading 
terms in l/Tgs, but rather follows from the behavior of 
the densities in the vicinity of the bubble wall. Thus, we 
conclude that Yb is dominated by the leading order con¬ 
tribution in l/Tss that results from retaining the finite- 
temperature masses in computing the statistical factors 
kj. 


F. Profile functions 


The CPV source term in Eq. (41) is proportional to 


the combination of VEVs which is convenient to 

rewrite in the following way: 


dll = CL22 = '^7 


<^12 = ^21 = 


Dh + {Dh + Dq){cq + 2ct) 
' D 

{cq + 2ct){Dh - Dq) 

' D 


rii = 


r2i = - 


[Dh + {cq + 2ct)Dq\{T M- + r^) 

knD 

{cq + 2ct)Dq{Tm- + ^h) 


knD 


ri2 — r22 — 0 


^1 = 


^2 = - 


Dh + {cq + 2ct)Dq 
D 

(cq + 2ct)Dq 




D 


H ’ 


D = Dh + 2(cq + 2ct)Dq . 


(53) 


The left handed quark density is found via ul = 4T + 5(5 
EQ). The electroweak sphalerons transfer this density 
into the net baryon asymmetry according to m 


vs{x)vj:{x) [ f 5 ( a ;) vE ( a ;) - 
= ^ [vE5(a;)]^ sin 2^Es(a;) pj:s{x), (55) 


where we have defined 




\JvUx) + vl{x), 


ks{x) = arctan 


vs(x) 


.(56) 


The detailed shapes of the profiles or equivalently 

across the bubble are unknown and their de¬ 
tailed calculation is beyond the scope of this paper. Anal¬ 
ogous calculation of the profiles in the MSSM have been 
performed in Ref. [53] and for the complex singlet ex¬ 
tension in Ref. [23]. In the MSSM case the role of two 
VEVs are playing Vu{x),Vd{x). Eor simplicity we adopt 
the shape of the profiles from MSSM calculations, but 
stress that further study is required for a more complete 
treatment. Thus, we use the following analytical form 
[4TH43l[54] of the profile functions 


Ub 


r 

—3—^ / dznL{z) ex.p 
Xw J — oo 


Vw J 


(54) 


Our expressions in Eq. ( [5T| are equivalent to the ri 
term in an analogous approximate solution in the MSSM 
case as given in Eq. (84) in Ref. [41]. Note that if one 
uses the zero mass limit values of kq = 2/ct = 2/cb, 
from Eq. (54), Eq. ([^, Eq. ([^ we see that the net 


baryon asymmetry equals to zero, as expected based on 


t'E5(^) = 24°l(r) 


^2az\ 


/^Es(x) = /3o(T) - -A/3 


1 + tanh 

1 — tanh 


/2az\ 


(57) 


with a = 3/2. Note that we have assumed the same wall 
thickness for the triplet and singlet vevs, an assumption 
that parallels the treatment of the doublet and singlet 
vevs in Ref. [23] . 
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Similarly to the MSSM we will assume additional sup¬ 
pression of the BAU due to small numerical value of 
A/3 0.015 [53] . Unlike the MSSM case our source term 

is proportional not only to but also to sin2/d(x), 

which in principle could lead to additional suppression of 
the total generated BAU. For concreteness we choose the 
value for Po{T) in the early universe to be 7r/4, so that 
sin2;do(T) = 1. Our results can be appropriate re-scaled 
after a comprehensive study of the bubble profiles has 
been completed. 


VI. RESULTS 

In this section we present the numerical results for 
bounds from existing and sensitivity regions from fu¬ 
ture generation EDM measurements, and also constraints 
from Higgs to diphoton decay and BAU. The question we 
are going to ask is whether it is possible to generate the 
observed BAU during the first step of the 2SPT described 
in this paper, in a way that is consistent with the exper¬ 
imental observations. 


A. EDMs 

From our previous analysis the dependence of electron 
and neutron EDMs on the parameters of the theory has 
the following form 

de = sin tan/3F(m/i3, m^), (58) 

dn = sin dx: tan /3 + CJ cot 13) F{mh ^, m^), 

where Cf, C 2 depend on the quark charges and nucleon 
matrix elements of up and down quark EDM operators. 
A similar expression to that for applies to the proton 
EDM with appropriate replacements C 12 ^ 12 - 

In Ref. [31] the nucleon tensor charges have been eval¬ 
uated on the lattice 

= 0.774(66), =-0.233(28). (59) 

Using the above results we obtain from Eq. ( [3^ the fol¬ 
lowing numerical values for the nucleon matrix elements 
of the quark EDM operators 

C“ = (3.5 ± 1.0) X 10“22cm, 

= (-24.2 ± 2.9) X 10“^2 
C“ = (-11.6 ± 3.2) X 10-22 cm, 

Cp = (7.3 ±2.1) X 10-22 cm. (60) 

Results based on sum rules and quark model can be 
found in the review m and are approximately factor 
of 2 smaller with larger error bars. We have also used 
the PDG [55] quark mass values rriu = 2.3 ± 0.6 MeV 
and rrid = 4.8 ± 0.4 MeV [64] and combined the rela¬ 
tive uncertainties from quark masses and tensor charges 


in quadrature, assuming they are uncorrelated. The re¬ 
sulting uncertainties in Eq. (60) range from 12% to 29%. 
Note, however that we have neglected possible contribu¬ 
tions from the strange quarks, whose contribution to dn 
may be as large as 35% with an uncertainty of similar 
magnitude}^. Consequently, to be conservative, we will 
take the uncertainty on to be 29%. 

The constraints on the CPV phase are shown in 
Figure as a function of tan^d,• The gray 
bands correspond to the current ACME electron EDM 
bound |de| < 8.7 x 10“^^ecm at the 90% confidence 
level. The light-red bands correspond to the projected 
neutron EDM sensitivity jd^l = 2.9 x 10“^^ecm, which 
is a factor of 100 times more sensitive than the cur¬ 
rent neutron EDM bound. The dark yellow bands corre¬ 
spond to a possible proton EDM search with sensitivity 
of \dp\ = 2.0 X 10“^^ecm. 

As expected for the first panel in the Figurethe elec¬ 
tron EDM bound on sin is a hyperbola as a function of 
tan /3, the neutron EDM bound becomes stronger at both 
large and small values of tan /3, consistent with paramet¬ 
ric dependence in Eq. (58). The second and third panels 
reveal the dependence of function F{mh ^, m^) in Eq. (58) 
on its arguments. Also as expected, the bounds become 
weaker as one increases the mass of either scalar that 
enters the two-loop contributions. In all three plots we 
indicate in the top left the numerical values used for the 
remaining two parameters, besides those that are varied 
in the plot. 

From the left panel, we see that the existing electron 
EDM bound yields tight constraints on the angle ex¬ 
cept at small values of tan p . A future EDM experiment 
with ^ 10“^^e cm sensitivity would cover this low-tan/3 
region. 

In addition we stress that the EDM bounds do not 
constrain the second CPV phase d^-, see Eq. ( [To] ) . We will 
exploit this feature below when we discuss the numerical 
results for constraints coming from the BAU. 


B. Higgs to diphoton 

Using the current LHC data from run I, the Higgs to 
diphoton signal is consistent with SM model and equals 

[321133] 

/^CMS = ( 61 ) 

Matlas = (62) 

These agree quite well and in our numerical analysis we 
use the combined result [56] 

MLHC = 1.15tS5;225 • (63) 

Note that the diphoton signal strength above corresponds 
to the production mechanism via gluon gluon fusion. 
Adding the corresponding vector boson fusion measure- 
ment 1.17t)!j| [56] will not change this result signifi¬ 
cantly because the error bar for it is a factor of 2 bigger. 
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Electron, Neutron and Proton EDMs 



Electron, Neutron and Proton EDMs 



Electron, Neutron and Proton EDMs 



m/j3(GeV) 


m2(GeV) 


FIG. 6: Electron EDM bound \de\ < 8.7 x 10“^^ecm (gray band represents the excluded region), projected neutron sensitivity 
|dn| = 2.9 X 10“^®ecm (pink band) and a possible proton EDM with sensitivity |dp| = 2.0 x 10“^®ecm (olive band). Widths 
of the pink and olive bands correspond to uncertainties in Eqs. (60). See text for more details. 


We see from equations Eq. (31)-Eq. (34) that all 
parameter-dependence of the Higss diphoton rate can 
be absorbed into where the gj: dependence on 

tanRe a 2 E is given in Eq. . The experimental 
constraints on the parameters tan/5,ms from the LHC 
diphoton data are shown in Figure We have chosen 
Rea 2 E = 1-07 — a value that is used in our BAU compu¬ 
tation below. The light-blue shaded region corresponds 
to the parameter space of the theory consistent with the 
LHC diphoton data at 90% CL. 


From Figure we find that for > 119 GeV any 
value of tdiiiP would in no tension with current LHC 
diphoton data. Thus we conclude that the run I LHC 
data on the Higgs diphoton decay is only weakly con¬ 
straining the parameters of our model; however we should 
keep in mind that run H data will place more stringent 
bounds. 


LHC bound from jiyy measurement 



FIG. 7: Light-blue region represents the parameter space of 
the theory that at the 90% confidence level is consistent with 
LHG diphoton data. 


C. Generation of the baryon asymmetry 

In this subsection we will show that for illustrative 
vev profiles vs{x),v^{x) adopted from MSSM studies 
and presented in Eq. ( [^ , and for typical parameters 
for the bubble wall width {L^j = 0.25/T) and wall ve¬ 
locity {vw = 0.05), it is possible to generate the ob¬ 
served BAU during the first step of the 2SPT. For 
the critical temperature and the values of the vevs of 
the triplet and the singlet in the early universe we use 
T = 123GeV, ue(T) = vs{T) = 76.3 GeV, motivated 
by benchmark studies of [12]. For diffusion constants we 
use Dq = 6 jT^ Dh = 110/T [49]. Again, a more compre¬ 
hensive study would require an explicit analysis of a more 
appropriate choice of profiles and bubble wall parameters 
needed for the scenario at hand. 

The BAU has a known resonant behavior [41] , reflected 
by a peak around some point in the parameter space. To 
see why this is the case, consider the CPV source term 

5'^, which for our problem is shown in Eq. ( [dlj ) . The 
function A to which the CPV source is proportional to 
is given in Eq. ( [4^ , and plotted in FigureThe reso¬ 
nance leading to the maximum magnitude for the CPV 
source corresponds degeneracy of the thermal masses of 
the two Higgs bosons m//^(T) = mH 2 {T) • It should be 
noted that similar resonant behavior is valid for the re¬ 
laxation rates Fm in Eq. ( [4Q| ), which to a certain degree 
mitigates resonance in the source [41]. The difference of 
mass squares of two Higgs bosons at finite temperature 
equals 


m 


Hi 




vl - Vt - ^ /3 - cot^ /?) 


rp2 

T- 


The BAU is, thus, maximal when tan (3 and mh 2 are cho¬ 
sen such that Eq. (64) equals zero. Note that the Yukawa 


parameters yt, yh also depend on tan [3 as indicated in the 
text below equation Eq. (38). 


As noted above, even though there exists no tree-level 
mixing between the two doublets due to the Z 2 symmetry 
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Resonance behavior of the function A, 



FIG. 8: Resonant behavior of BAU dependence on the pa¬ 
rameters of the theory. See text for more details. 


of V (i^i, i^ 2 ) (^12 = 0), at the finite temperature, an off- 
diagonal mass term ^m^(T)i7|i72 + h.c. is generated with 

5 im?‘{T) = — h ^^24 • non-negligible, such term 

would induce complicated flavor oscillations [25], requir¬ 
ing proper treatment beyond the VIA. For consistency we 
estimate the error made by neglecting these flavor oscil¬ 
lations by defining parameter 0 r\j 

that characterizes the magnitude of flavor oscillations 
Pose ^ [Is]- Therefore, if we choose to maximize 

the BAU and tune the parameters to be at the resonant 
point, the approximation of neglecting flavor oscillations 
becomes arbitrarily bad, or ^ = oc. For practical pur¬ 
poses of selecting benchmark scenarios, one should select 
parameters away from the resonance, keeping the mag¬ 
nitude of 0 from oc at reasonably controlled value. For 
the benchmarks that we study below the value of is 
included into the Table in. 


With these considerations in mind, we select two 
benchmark scenarios A and B that yield the observed 
BAU while respecting the present EDM bounds and 
avoiding significant flavor oscillations. We set Ss = 0 
for benchmark A and = 0 for the benchmark B. The 
remaining parameters, along with the thermal masses, 
relaxation rates in the broken phase, and the values of 
the EDMs and BAU are summarized in the Table [ll| [65] . 
As one can see from that table, the neutron EDM is cur¬ 
rently consistent with benchmark A, however with a fac¬ 
tor of 100 increase in the sensitivity this benchmark will 
be probed and possibly ruled out. On the other hand 
benchmark present and future EDM searches have no 
sensitivity to the benchmark B. 


The profiles of the relaxation rates, sources and 
particle-antiparticle asymmetries for benchmark point A 
are illustrated in the Eigurej^ The first panel represents 
the dependence of two effective relaxation rates Fn and 
r 2 i entering the approximate equations in Eq. (52) on 
the co-moving distance from the bubble wall, 2 ; = 


TABLE II: Table of parameter values, thermal masses, relax¬ 
ation rates, values of EDMs and BAU for benchmarks A, B. 


Parameter 

A 

B 

mh 2 [GeV] 

180 

180 

nihs [GeV] 

650 

650 

tue [GeV] 

130 

130 

mH+ [GeV] 

300 

300 

tan {3 

0.4 

0.4 

sin (5 e 

0.38 

0 

sin 5 s 

0 

-0.4 

Re a 2 s 

1.07 

1.07 

\ci2s\ 

2.0 

2.0 

hA 

0.8 

0.8 

m,ST) [GeV] 

106 

106 

mt^{T) [GeV] 

133 

133 

[GeV] 

63 

63 

muAT) [GeV] 

82 

82 

niHAT) [GeV] 

215 

215 

mjr) [GeV] 

91 

91 

ms{T) [GeV] 

200 

200 

Tifi [GeV] 

5 

5 

Th, [GeV] 

5 

5 

Ty [GeV] 

1.9 

1.9 

Th [GeV] 

0.31 

0.26 

Tss [GeV] 

0.41 

0.41 

rM+ [GeV] 

-0.084 

-0.080 

Tm- [GeV] 

0.48 

0.46 

de/ (10 '^'^eem) 

5.8 

0 

dn/ (l0“^®ecm) 

-3.3 

0 

dp/ (l0“^®ecmj 

6.0 

0 

yapprox/ (10-1'^) 

9.3 

9.1 

Vyil 

8.6 

8.4 


0.0073 

0.0066 

M 77 

0.86 

0.86 


with 2 : < 0 corresponding to the symmetric phase (out¬ 
side of the bubble) and 2 : > 0 to the broken phase (inside 
the bubble). The second panel represents the two sources 
81,82 entering the approximate formula. Note the quali¬ 
tative features that the relaxation rates F are zero in the 
symmetric phase and look like a step function across the 
bubble wall, while the CPV sources are zero everywhere 
except within the bubble wall. The smallness of the mag¬ 
nitude of the 82 compared to 81 can be understood from 
explicit formulae for them in Eq. (53) and the fact that 
numerically Dq ^ Dh- 

The third panel represents all the four densities 
T,Q,H,h in colors red, green, blue black respectively. 
The thin solid lines correspond to the full numerical so¬ 
lution to the transport equations while the thick dashed 
lines to the approximate method. Einally, the last panel 
represents the left handed quark density = 4T + 5(3 
[50] , which is converted to the net baryon asymmetry via 
electroweak sphalerons m , see Eq. ( |5^ . Eigurej^ serves 
as an illustration how particle-antiparticle densities are 
distributed inside and outside of the bubble in the early 
Universe. The last panel also shows that the approximate 
method to solve the transport equations is in a reason- 
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-0.006 -0.004 -0.002 0.000 0.002 0.004 0.006 


Model A: profiles of S\,S2 



z{GeV-’) 



z(GeV-’) 

Model A: profiles of 



FIG. 9: Profile functions and the solution to the boundary problem for benchmark A. See text for more details. 


able agreement with the full numerical solution, except 
in the vicinity of the bubble wall. While the agreement 
between the full and approximate solutions in the third 
panel might lead one to expect similar agreement in the 
last panel, we note that the LH density = 4T + 5(3 is 
subject to large cancellations. Consequently, any small 
differences between the full and approximate solutions in 
the third panel become magnified in ni, (note the order 
of magnitude smaller vertical scale in the fourth panel as 
well). Thus, it is not surprising that expectations based 
on the approximate solution are fairly well reflected in 
the final results for Yg. 

The resulting baryon asymmetry as a function of the 
CPV phases and Ss and 171^2 are shown in Fig.[^ In 
the left (right) panel we set ^5 = 0 (^e = 0). We also 
indicate the present EDM constraints and prospective fu¬ 
ture sensitivities, which appear only in the left panel since 
for = 0 the interactions in V (i^i, 7^2, E, S') generate no 
elementary fermion EDMs through two-loop order. The 
present electron EDM bound \de\ < 8.7 x cm 

obtained by the ACME collaboration [57] excludes the 
shaded region above the horizontal black line. The hor¬ 
izontal pink band indicates the reach of the future neu¬ 
tron EDM search underway at the Eundamental Neu¬ 
tron Physics Beamline at the Spallation Neutron Source 
that has a goal sensitivity of \dn\ = 2.9 x 10“^^ecm. 
A possible future proton EDM search with a sensitivity 


of \dp\ ^ 10“^^ecm would cover the entire BAU-viable 
region of the left panel. Eor illustrative purposes, we 
show the reach with a proton EDM experiment having 
\dp\ = 2.0 X 10“^^ecm sensitivity with the olive band. 
The widths of the proton and neutron EDM bands cor¬ 
respond to 29% error, as in Eig. The green bands cor¬ 
respond to the parameters for which the observed baryon 
asymmetry is generated. The discontinuity in the slope 
at mh 2 ^ 350 GeV results from crossing a kinematic 
threshold in the three-body, particle number changing 
rates. 

These results indicate that it is possible for the ob¬ 
served BAU to be generated during the first step of the 
two-step EWPT. The present electron EDM bound ex¬ 
cludes a portion of the BAU-viable parameter space asso¬ 
ciated with the CPV phase while future nucleon EDM 
searches could probe most or even all of this sector of the 
model. On the other hand, the source of the BAU asso¬ 
ciated with the singlet-Higgs operator a 2 sHlH 2 S‘^ + h.c. 
is immune from these present and future EDM probes. 

It is also interesting to ask how these statements vary 
with the other parameters in the theory, particularly 
tan/3. To that end, we show in Eig. the tan/3- 
dependence of the CPV phases and EDM sensitivities. 
We restrict our consideration to tan/3 ^ 5, as for val¬ 
ues above this region, one must include explicitly the ef¬ 
fects of bottom quark Yukawa rates that are enhanced as 
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Model A: 


Model B: ^2=0 




FIG. 10: Constraints on the CPV phases fe (left panel) and Ss as a function of the mass mh 2 - Solid green band is 
consistent with the observed BAU: Yb = (8.59 zb 0.22) x 10“^^. For Ss = 0 (left panel) the current electron EDM bound 
\de\ < 8.7 X 10~^^e cm excludes the shaded region above the horizontal black line. Sensitivities of a future neutron EDM 
\dn\ < 2.9 X cm[53 and possible proton EDM with \dp\ < 2.9 x 10~^^e cm are indicated, respectively, by the pink and 

olive shaded horizontal bands. Widths of the bands correspond to uncertainties in Eqs. ( |60[ ). See text for more details. Present 
and future EDM searches have no sensitivity when = 0 (right panel). See Table [n| parameters corresponding to benchmarks 
A (magenta circle) and B (blue square). Eor each panel, we fix these values and vary only the two parameters shown. 


tan^/5[47] in the type II 2 HDM. Note that for Model A, 
the present electron EDM bound restricts one to values of 
tan/3 below unity. For Model B, where the EDM places 
no constraint on the CPV phase 6 s ^ we observe that 
consistency with the observed BAU requires tan^d < 3. 
Based on experience with the MSSM, it is likely that 
in the large tan (3 regime for both models there may be 
significant cancellations between effects associated with 
top and bottom Yukawa rates, leading to an even smaller 
Yb and the requirement of larger CPV phases. Conse¬ 
quently, we expect that this regime will not be viable, 
and we defer a detailed study of this regime, as well as 
an analysis of the type I 2HDM realization of our model, 
to future work. 

In principle, it would also be interesting to explore the 
mE-dependence. In the present set-up, larger values of 
ms would not be consistent with the two-step EWSB sce¬ 
nario according to the analysis of Ref. [ 12 ] . For smaller 
values of ms, the EDMs induced by (3s become larger, 
leading to ever more severe constraints for Model A. We 
observe that allowing (a) the coefficient of the op¬ 

erator to be non-vanishing and (b) (S) 7 ^ 0 at T = 0 
would open the possibility of larger values of the T = 0 
triplet mass, thereby in principle weakening the EDM 
constraints on Model A, while preserving the viability of 
the two-step EWSB transition at finite T. We again defer 
a detailed study of this possibility, including the impact 
on asymmetry generation, to future work. 


VII. CONCLUSIONS 

The possibility that EWSB occurred in multiple steps 
opens a new pathway for weak scale baryogenesis. In 


the two-step EWBG paradigm, the BAU is generated 
by CPV dynamics during a first order transition to an 
EWSB-vacuum that precedes the final transition to the 
Higgs phase of today’s universe. As outlined in Section 
|n| there exist multiple possibilities for the CPV inter¬ 
actions that could make this new pathway effective. In 
this study, following the previous work in Ref. [ 12 ], we 
have illustrated one subset of these possibilities, focusing 
on renormalizable interactions in the scalar sector that 
directly generate CPV-asymmetries in the SM Higgs sec¬ 
tor via Higgs portal interactions. In this example, BAU 
generation requires the presence of two fields that ob¬ 
tain space-time varying vevs during the first EWSB step 
and two fields in the Higgs sector that mix due to CPV- 
interactions with one or both of these vevs. For the for¬ 
mer, we choose a real triplet Y, whose vev breaks EW 
symmetry during the first step, and a real singlet. For 
the latter, we employ a type H 2HDM, wherein neither 
of the neutral doublet fields obtain vevs during the first 
step but do so in the second step (see Fig. Q. The real 
singlet provides an example of a “partially excluded” sec¬ 
tor that interacts with the SM solely via the new CPV 
interaction. 

Several generalizable features emerge from our illustra¬ 
tive model study. Most significantly, the two-step EWGB 
paradigm appears to be a viable mechanism for creating 
the BAU and does not appear to require fine tuning of 
parameters. Moreover, while a portion of the CPV dy¬ 
namics is accessible to present and future EDM searches, 
those associated with the partially secluded sector are 
not. Direct searches for the new scalar states could ei¬ 
ther discover exclude the ingredients necessary for this 
scenario, but direct tests of the CPV interactions are lim¬ 
ited to those involving non-singlet fields, at least for the 










Model A: 6s=0 


Model B: (^2=0 
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tan p tan p 

FIG. 11: Constraints on the CPV phases (left panel) and 6s as a function of the tan/3 for hxed mh 2 = 180 GeV. Various 
regions and curves have the same meaning as in Fig. |10[ 


foreseeable future. On an experimentally more positive 
note, should the CPV responsible for the BAU involve the 
non-secluded sector, one could anticipate non-vanishing 
signals in future EDM searches. 

Looking ahead, it would be interesting to explore both 
other specific realizations of two-step EWBG as well as 
to study the present example with greater comprehen¬ 
siveness. The latter analysis would include considera¬ 
tion of the type I 2HDM; allowing for a non-vanishing 
coupling that could yield larger triplet masses con¬ 
sistent with the strong first order EWPT during the first 
step; carrying out a detailed study of the bubble profiles 
and wall velocities; and ultimately going beyond the VIA 
along the lines of Refs [25l [26] . 
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TABLE III: Thermal masses. 


Thermal mass 6m%yi/T‘^ 
q93 + ^92 + iQy^ 


Ql 

tn 

6r 

Hi 

H2 

E 

S' 


PI + R 
hi + ^9l + H 
+ ^9l + hi + T 


_ 1 _ As _i_ A 4 

“C 6 “C 12 


Appendix B: Thermal masses 

Thermal masses for the fields ^ 1 , ^ 2 , 5" 

are summarized in the Table mi 


Appendix A: EDM loop integrals 


Here we summarize the loop functions needed for the 
EDM calculation. 


Appendix C: Analytical integration for the 
boundary problem 

Solving the boundary problem for the densities vary¬ 
ing across the bubble wall is equivalent to considering N 
coupled linear second-order differential equations 



f{x, y) 
9{x,y) 


vf{x) 

y-x 

ygjx) 

y-x 


[1 — 2x{l — x)] In 
x{l — x) — z 

x(l-x) 

_^_ 

x{l — x) — 

xfjy) 

x-y' 

xgjy) 

x-y' 


(Al) 


y-(z) + tti y[{z) + bij{z)yj{z) = Si{z), (Cl) 

where index i = 1,..., A is fixed and not summed 
over. Eor our problem ^ are constants while 

hij{z) ^ Vij{z)lDi vary across the bubble wall and dif¬ 
fer in the symmetric and the broken phases. However 
away from the bubble wall | 2 :| ^ Lyj the functions bij{z) 
converge to constant numbers which we define as: for 
z > 0, bij{z) ^ 6F and for z < 0, bij{z) ^ b\j , where su¬ 
perscripts r, I stand for “right” and “left” corresponding 
to the “broken” and “symmetric” phases respectively. 
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It is c onvenient to reduce the system of equations 
Eq. (Cl) to the 2N x 2N first order coupled differential 
equations by introducing a new variable 


Yiz) = 


- ( Viiz) 


y'iiz) 


(C2) 


which is a column vector of the size 2N. In the column 


notation the set Eq. (Cl) reduces to 


Y\z) =A^Y{z)^S{z). 


(C3) 


In the equation above ^4 is a 2N x 2N dimensional matrix 
and S' is a 2N dimensional column vector which in the 
block-diagonal form equal to 

where a and b are both N x N dimensional matrices with 
a = diaga^ and b matrix has matrix elements equal to 
bij and we emphasized again that we are working un¬ 
der the assumption of no z dependence of the matrix 
A. A/"—dimensional column vector s{z) has elements 

equal to Si{z). Define U to be a matrix that has its 
columns consisting of the eigenvectors of the matrix A 
and assume that it diagonalizes the matrix A according 
to U~^AU = ^diag = diagA/c, where Xk are the eigenval¬ 
ues of the matrix A. Because we assume that A has no 
dependence, likewise we obtain that t/, Xk have no 
dependence. In this case the transformation Y = U~^Y 
leads to a simple uncoupled set of 2N linear differential 
equations of the first order 

n = XkYk + Sk{z), S = U-^S{z). (C5) 


The solution in terms of the initial conditions for Y is 


Yk{z)=e^^^ 



(C6) 


Note that in the symmetric phase {z < 0) and in the 
broken phase {z > 0), the matrices U^Xk/S are differ¬ 
ent. We will assume a superscript l^r where appro¬ 
priate to identify on which side of the bubble wall we 
are studying the solution to the boundary problem. The 
boundary problem requires Y (—oc) = Y (+oo) = 0 which 
easily translates into the boundary condition for the ro¬ 
tated variables T(—oc) = T(+oo) = 0. Note that while 
the formulation of the boundary problem requires that 
the column vector Y(z) is continuous across the bubble 
wall T(—0) = T(-FO), the same is not true for Y{z): 
T(—0) 7 ^ T(+0). The necessary conditions for the 

boundary problem to have a solution is 


/ dte S'fe(i), 
Jo 


U(+0) = - ‘ 

p — OO 

Yk{-0) = - dte-^'Y 
Jo 


lim 

z^-\-oo 

lim 

z^ — oo 


[ dte^’^^^-*'> Skit) =0, 
Jo 

r Skit) =0, 

Jo 


\l > 0, (C7) 
Ai <0,(C8) 
A^<0, (C9) 
Ai>0.(C10) 


Note that all integrals in equations Eq. ( |C7[ )-Eq. ( |CIQ ) 
are finite which follows from integrability of Sk{t) near 
oc. Simultaneous solution of the equations above to¬ 
gether with continuity condition at ^ = 0 solves the 
boundary problem. If such solution does not exist then 
the boundary problem has no solution. 

In practical applications it is often used the following 
VEV profile functions 


viz) = 


- tanh ^ 


P'iz) = 


2Lw cosh ^ 


{C^ii) 


with Si ^ v{z)^fi'{zY • Eor example in the SUSY ex¬ 
ample a = 2^ = 1. In the present paper we have 
(a = 4, /3 = I. The following master integral we eval¬ 
uate analytically 


=/■ 


7(z,A;a,/3)= / dte’"* [mV 


- 0 


2ia + l3)-L^X 


(C12) 


Lw X 


X 2 E 1 I ot 2/3 ,0 4-/3— —-— , I 4“ o 4“ /d — 


LwX ^ 


Consistently with our assumption that the relaxation 
rates are approximately given by a step function across 
the bubble wall it is safe to assume that |z| in the 

equation above. Indeed the BAU is generated far away 
from the bubble wall and this is an extremely reliable 
approximation. Thus by taking limits of the hypergeo¬ 
metric function at infinity obtain 


/(-hoc. A; o,/3) 

2Ulr^ 


2(o; -\- d) — LwX 


r (1 + a + /3 - r (^ + 
r (a + 213) 


- 2 F 1 (a + 2/3,a + ^-^,l+a + /3-^;-l 


r ^ 2/3 

for A > — . 

-L^w 

In the opposite limit ^ —>■ —oo we find 
/(-oo, A;a,/3) 

^ 2Ul-^ 

2(0 + /3) — LwX 


(C13) 


(C14) 


X 2 F 1 (a + 2/3,a + /3-^,l + a + /3-^;-l 


for A < 


2(0 + /3) 


The equations Eq. ( CI3] ) and Eq. (CI4) are directly ap¬ 
plicable in Eq. (C7) and Eq. (C8). Einally we also es¬ 
tablish the following limits 

lim e^^/(z. A; o,/3) = 0, for A < 0, (CI5) 

2:^-|-oo 

lim e^^/(z. A; (a,/3) = 0, forA>0.(CI6) 


T he eq uations above demonstrate that Eq. (C9) and 
Eq. (CIO) automatically are satisfied. 
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